Spin-gap effect on resistivity in the t-J model 
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We calculate the spin-gap effect on dc resistivity in the t-J model of high-T c cuprates by using 
the Ginzburg-Landau theory coupled with a gauge field as its effective field theory to get p(T) oc 
T{1 — c (T* — T) d }, where T* is the spin-gap onset temperature. By taking the compactness of 
massive gauge field into account, the exponent d deviates from its mean-field value 1/2 and becomes 
a nonuniversal T-dependent quantity, which improves the correspondence with the experiments. 
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The idea of charge-spin separation (CSS) by Anderson 
[jjj which accounts for the anomalous behavior of various 
normal-state properties of high-T c cuprates M allows us 
to treat holons and spinons introduced in the slave-boson 
(SB) mean-field theory (MFT) of the t-J model as quasi- 
free particles. Fluctuations around MFT are described 
by gauge fields coupled to holons and spinons, the effects 
of which may be calculated in perturbation theory. Ac- 
tually, by using a massless gauge field, Nagaosa and Lee 
H obtained the dc resistivity p(T) oc T, which agrees 
with the experiment for T > T* m where T* is the onset 
temperature of spin gap. 

For T < T* , the experimentally observed p(T) reduces 
from this T-linear behavior |||(| . Because the gauge field 
is expected to acquire a mass tua in the spin-gap state, 
this reduction could be understood as a mass effect; the 
fluctuations of the gauge field become weaker and the 
scatterings between holons (the carriers of charge) and 
gauge bosons are reduced. Actually, in the previous pa- 
per , we obtained the following result M ; 



P(T) 
X(T) 



2e 2 ns 



1 + X{T) - y/{l + X(T)} 2 -1 



1 



87rns(T)' m 



1 

rriF 



2m,B 



(1) 



where rriFiB) is the spinon (holon) mass, /s(— Hb) = 
[exp(— Ppb) — 1] j with the holon chemical potential 
fJ-B, ns is the holon density, and Ub = Ub / /b(— A*b)- 
The factor in the square brackets represents the reduction 
rate due to m\(T). By assuming the behavior m^(T) oc 
(T* - T) d {d > 0) near T*, and ignoring the weak T- 
dependence in m(T) and ub{T), we get 

p(T)<xT[l-c(T* -T) d }. (2) 

The MF value d = 1/2 is excluded since the experiment 
shows smooth deviations from the T-linear form, which 
requires d > 1. In this letter, we calculate uia by setting 



up the effective field theory and taking the compactness 
of gauge field into account, finding that d is not a univer- 
sal constant but has the T-dependence. 

The effective theory is the Ginzburg-Landau (GL) the- 
ory L(X,Ai) coupled with a gauge field Ai, where the 
complex scalar A represents the e?-wave spinon pairing. 
Halperin, Lubensky, and Ma j^] considered a similar sys- 
tem in 3D, where Ai corresponds to the electromagnetic 
field and A to the Cooper-pair field. They calculated 
the effect of Ai on A, which converted the second-order 
phase transition to first-order. More recently, Ubbens 
and Lee jl(J calculated the one-loop effect of Ai in the 
SB MFT of the 2D t-J model, and concluded again that 
the pairing transition at T* becomes first order. How- 
ever, their T* appears below the superconducting tran- 
sition temperature T c , so they concluded that the spin- 
gap phase is completely destroyed by gauge-field fluctu- 
ations. In the present study, we take the compactness 
of Ai into account, which originates from the t-J model 
defined on the lattice itself and gives rise to interactions 
like A 2 cos(2a^4i), where a is the lattice constant. Even 
in the CSS state, it generates nontrivial vertices that are 
missing in the usual treatments (!(]] which use X 2 Af. We 
find that the periodic interaction stabilizes the system so 
as to have a phase transition above T c . 

For the 2D GL theory coupled with a gauge field, 
Nagaosa and Lee 11 argued that vortices put a phase 
transition into a crossover. On the other hand, the 3D 
system has a genuine phase transition. This is sup- 
ported by Monte Carlo simulations and other studies jl2) . 
In the low-temperature phase (spin-gap phase, Higgs 
phase), vortex loops do not proliferate, while in the high- 
temperature phase, they do. The former phase is well 
described by the usual order parameter A, while the lat- 
ter is described by the disorder parameter that measures 
the vortex-loop density. We assume a small but finite 
three-dimensionality with anisotropy a > (a = 1 for 
3D and a — for 2D), so there takes place a genuine 
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phase transition at T* in the present model. 

For a sufficiently small a, the 3D critical behavior is 
observed only in the small interval in T near T* which 
vanishes as a — > 0. Beyond this, calculations in the pure 
2D system (a = 0) should give a reliable result according 
to the general theory of critical phenomena. There is a 
good example of this; the antifcrromagnctic transition of 
cuprates at T = Taf for which the 3D behavior 

takes place for \T ~ T AF \ < 0(1/ 1 lna| 2 ). Beyond this 
interval, the 2D results fit the experimental data well. 

Let us start with the SB t-J Hamiltonian given by 
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where /„ is the fermionic spinon operator with spin a 
(=T,I) at site i of a 2D lattice, and b x is the bosonic 
holon operator. The direction index i (— 1,2) is also 
used as unit vectors. p>B,F are the chemical potentials to 
enforce {b\.b x ) = 5j_Z) ct (/J ./ :c <t) =1-8 where S is the 
doping parameter fll4j| . We introduce the complex auxil- 
iary fields Xxi and X X i on the link (x, x + i) to decouple 
both t and J terms as 

x, i 

-{x*i(ljY,f*+i°f*°+ tb Ui b *) + H - c -} ( 4 ) 

Xxi describes hoppings of holons and spinons, while X X i 
describes the resonating-valence-bond (singlet spin-pair) 
amplitude. We shall treat their radial parts as MF's. 

In path integral formalism, the partition function 
Z(f3) [/? ee (kvT)- 1 } is given by 



Z= I {db}[df][d X ][dX\ exp(-S), 
/o 



fxafxaj + HmF 



(5) 



where r is the imaginary time and b x = db x {r) j 3t , etc. 
Let us consider the low-energy effective theory at tem- 
peratures T <C Tcss, where less is the critical tempera- 
ture below which the CSS takes place as a deconfincmcnt 
phenomenon |fl5|| . It is reasonable to translate the lat- 
tice variables to the continuum fields like f xa — > af a (x), 
etc. using the lattice constant. The Hamiltonian of the 
continuum field theory is given by 
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where plb,f are the chemical potentials for the continuum 
theory, f a (k) is the Fourier transform of f a {x), and 
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SG (k, q) = 2(1 - 8)(^-^)x d {q) - 25 X s (q). (7) 
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To obtain H, we modified the dispersions of holons and 
spinons from cosine form to the quadratic one. Di = 
di — iAi is the covariant derivative with the gauge field 
Ai{x). Here we introduced Ai(x) by the correspondence 
Xxi X e:>0 P[i a Ai(x)], where x is the radial part of Xxi- 
We ignore the fluctuations of \xxi\ since T <C T x , where 
T x is the onset temperature of x- I n the SB MFT, x is 
estimated at small 6's as 
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if the spinon pairing X X i is neglected. 

To obtain the effective action of Xi and Ai, b and Jo- 
are integrated by the standard bilinear integrations. This 
procedure generates dissipative terms of Ai and Ai. The 
most singular contributions to Z from the integrations 
over A; and Ai come from their static (r-independent) 
modes, so we keep only the static modes in the effective 
Lagrangian density which is given up to the fourth-order 
in fields and derivatives by 

L cB = a s \X s \ 2 + a d \X d \ 2 + AbS 4 \X s \ 4 + h (1 - «5) 4 |A d | 4 
+ 2b5 2 (l-S) 2 (A\X s \ 2 \X d \ 2 + X 2 X 2 + X 2 d Xl) 
+ c ^(25 2 |P 4 A S | 2 + (1 - S) 2 \V t X d f 



c (5(1 - S) (ViXsVxXd - V 2 X s V 2 X d + H.c. 
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where m is defined in Eq. ([!]), and 
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where 7 is the Euler number, is the cutoff of the 
spinon energy [£ = fc 2 /(2mj^) — /if] in the one-loop in- 
tegrals representing spinon pairings, and is estimated as 
uj\ ~ 0(/ii?). From the potential energy of Aj(a;), the sys- 
tem favors the d-wave state at small 8's, and the s-wave 
state at large <5's. Let us focus on small S's by parameter- 
izing \%(x) = Xexp[i8(x)], X%(x) = —Xexp[i9(x)]. Here 
we introduced A, the spin-gap amplitude, for the radial 
parts of Xi(x), ignoring their fluctuations. Then we have 
the effective Lagrangian density, 



L c ji — L\ + La, 
Lx = a d X 2 + \b (1 - 8fX\ 



(11) 
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From the above La> the MF result is obtained as 
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where T\ is the critical temperature below which A de- 
velops, and should be identified as T*. The second result 
is reliable at T near T\. 

Let us consider the effect of La on L\ by integrating 
over Ai . We have treated Ai as a noncompact gauge field 
although it was originally compact. This procedure is ap- 
propriate for the kinetic term of Ai, because we consider 
the region T <C less- We will respect the compactness of 
Ai and the angle-nature of by considering the following 
new Lagrangian Lb with the periodic mass term; 
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where we introduced the Proca (massive vector) field 
Bi = Ai - di6/2. (Fij = d z Bj - djBi.) Let us 
take the unitary gauge. Then the integrals reduce as 
[dd][dAi\ = [d9][dBi] [dBi]. 

Let us estimate the gauge-field mass by the variational 
method. We choose the variational Lagrangian L' B for 
Lb as 
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where is a variational parameter. The variational 
free energy density Fb = F' B + (Lb — L' B )'is given by 
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where g c , the momentum cutoff of Bi, is 0(a _1 ). Wc 
have omitted the higher-order terms of 0(m A A / q^). Note 
that we took the definition of the propagator at the same 
point and the trace of a functional operator O as 



{B i {x)B j {x)) = UmJB^B^y)), 
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By minimizing Fb{ttia), we get 
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We note that the fluctuations of Ai do not affect the MF 
result (|l^) as long as d(T) > 1 since then the order of 
the corrections becomes higher than A . Thus the gauge- 
field mass m A (T) starts to develop continuously at T\ as 
m A (T) oc (T x ~T) diTx K That is, the exponent d is neither 
1/2 nor a constant, and drastically changes especially 
when T\ ~ Ta, where Ta is a root of the equation Ta = 
Q(Ta), at which d(T) diverges. This is in strong contrast 
with the noncompact case [|10| . 

If we write q 2 — ea~ 2 , we have m 2 A (T) — q 2 z (1 — 
T/T x ) 2d{T) with z = [16vr(l - 8)fh/ (em F )] 2d(T) ■ A 
straightforward estimation, irq 2 = (2ir/a) 2 , by keeping 
the area of momentum space, gives e — 4m. However, 
this gives rise to a nonrealistic curve of p(T) that devi- 
ates from the T-linear behavior and decreases too rapidly, 
due to the large factor z ~ 4 2d ( T ). Thus we regard q c to 
be a parameter of the effective theory, and choose e so 
as to obtain a reasonable p(T). For example, we require 
z = 1, which implies e = 167r(l — 8)fh/rriF- 

Finally, we need to consider the renormalization effect 
of the hopping parameter t. We assume that the 3D 
system exhibits Bose condensation at the temperature 
scale of Tb — 27m b I tub — AittxS, and regard Tb to be 
the observed T c in the lightly-doped region. Since t ~ 0.3 
eV gives rise to Tb ~ 3000 K at 8 ~ 0.15, one needs to 
use an effective t* ~ 0.01 eV in place of t so as to obtain 
a realistic T c ~ 100K Q. 

We show in Fig^ the phase diagram with the spin- 
gap on-set temperature T\, Ta at which the mass ex- 
ponent d diverges, and Tb- In Fig||, we plot p(T). As 
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explained, the curves reproduce the experimental data 
much better than those with d = 1/2 of the MF result, 
showing smooth departures from the T-linear curves, i.e., 
d(T\) > 1 for the region of interest, 0.05 < S < 0.15. 

The present results of p(T) support that our treatment 
of gauge-field fluctuations by the variational treatment of 
compactness is suitable to describe the spin-gap state in 
the t-J model, although more investigation is certainly 
necessary. 
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FIG. 1. Mean-field phase diagram of the t-J model. Tb 
is the Bose condensation temperature. T\ is the spin-gap 
onset temperature. Ta is the root of T = Q(T) at which 
d(T) diverges. We chose t* = 0.01 eV, J = 0.15 eV, and 
lu x = nJ X /(2e->). 
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FIG. 2. Resistivity p(T) in h/e 2 for several S's with the 
parameters chosen in Fig.l. The dotted lines represent the 
case of X(T) = in (1). The exponent d(T\) decreases as 
16.4, 4.8, 2.8, 2.0, as 5 increases. 
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